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FHEFACE 

nils  ^feaoraDdum  Is  Intended  to  be  a  fundaaental  old  In  studies  of  orblt> 
al  aedianles.  The  relatively  sliqpde,  straightforward  eoevutatloaal  procedures 
described  herein  can  be  used  to  establish  approxlaate  heliocentric  transfer 

I 

orbits  and  their  elenents.  On  the  basis  of  these  approadaate  orbits,  exact 
Interplanetary  orbits  can  be  calculated. 

These  cooputational  procedures  will  be  of  use  to  persons  engaged  In 
space>survelllance  activities  and  in  Interplanetary-mission  lOannlng. 


SUMMARY 
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This  Meinorandum  presents  and  discusses  some  methods  for  establishing 
heliocentric  interplanetary  transfer  orbits.  The  four  basic  methods  and 
tlieir  variations  can  be  used  to  establish  oi’bits  having  specified  transfer 
angles,  transfer  times,  hyperbolic  excess  velocities,  or  heliocentric  de¬ 
parture  velocities.  Each  method  consists  of  a  step-by-step  computation  pro¬ 
cedure  which  utilizes  the  equations  of  two-body  motion  and  appropriate  trigo- 
noEctric  relations  to  establish  the  desired  transfer  orbit. 

Each  of  the  methods  for  establishing  a  desired  transfer  orbit  requires 
an  Iterative  process.  Thus,  the  methods  are  best  applied  by  using  a  large- 
scale  digital  computer.  In  this  way  numerous  orbits  can  be  established  and 
the  orbit  which  is  optimum  for  some  requirement  can  be  selected.  None  of 
the  methods  permits  a  direct  analytical  determination  of  on  optimum  orbit. 
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LIGT  OF  SYMBOLS 

semimajor  axis 

semimajor  axis  and  eccentricity  of  EarUj's  orbit 
semimajor  axis  and  eccentricity  of  the  transfer  orbit 
eccentric  anomaly 

nth  positions  of  the  Earth,  vehicle,  and  destination 
planet 

eccentricity 

inclination  of  destination  planet  orbit  plane  to  the 
ecliptic 

inclinations  of  the  transfer  orbit  to  the  ecliptic  plane 

arbitraiy  constant 

gravitational  parjimeter 

mean  angular  rate  of  the  vehicle  in  orbit 

serailatus  rectum  of  destination  planet  orbit  and  Earth 

orbit 

distance  from  center  of  force  to  body 
distance  from  the  Sun  to  the  Efirth  at  departure  time 
distance  fixim  the  destination  to  the  Sun  at  orilval 
time 

perihelion  distances  of  the  orbits  of  the  destination 
planet  and  Earth 
an'ival  time 
departure  time 


transfer  time 


M-2881 

X 


V  ,  =  heliocentric  departure  velocity  of  the  vehicle 
sd 

=  hyperbolic  excess  velocity 
V  =  true  anomaly 

V  -  V  =  vehicle  transfer  angle 
s  sd 

v^dj  "  true  anomaly  of  the  vehicle  at  departure  and  arrival 
a  =  any^e  between  the  velocity  vectors  of  Earth  jmd  the 
vehicle  at  departure  time 

7^^  =  angle  betwen  the  heliocentric  velocity  vector  of  Earth 

and  norma],  to  the  Earth-3un  line,  measured  in  the  ecliptic 
plane 

7  ,  =  angle  between  the  vehicle's  velocity  vector  and  normal 
gcl 

to  the  ■i-'ehicle-.jun  line,  measured  in  the  transfer-orbit 
plane 

cL'p  =  argurent  of  perigee  of  the  destination  planet 
T]  =  mean  angular  rate 


=  iieliocentric  longitude  of  the  Earth  at  departure  time 

A  ,  X  =  heJ.iocentric  longitude  and  latitude  of  the  destination 
pa-^  pa 

planet  at  arrival  time 


A^^  =  heliocentric  longitude  of  the  vehicle  at  departure  time 
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I.  IWTROIXJCTION 

Interplanetazy  transfer  orbits  may  be  establl^ed  subject  to  one  or 

more  of  several  possible  constraints.  For  example,  one  may  vant  to  find 

all  orbits  vhlch  have  the  same  characteristic  depairture  velocity,  or  a  fixed 

(l  2) 

time  of  transfer  may  be  the  most  important  consideration.^  ’  Ihere  are 
several  other  considerations  such  as  arrival  date,  true  anomaly  of  the  de¬ 
parture  point,  etc. 

Ibis  memorandum  presents  and  discusses  some  coeqputatlon  procedures  for 
detenalnlne  the  elements  of  elliptical  Interplanetary  transfer  orbits  sub¬ 
ject  to  some  of  the  above  considerations.  Specifically,  several  different 
approaches  are  presented  and  the  cooputatlon  procedure  for  each  Is  given. 

For  all  the  cccputatlon  procediires  the  transfer  orbit  Is  assumed  to  extend 
between  two  massless  points  \dilch  coincide  with  the  center  of  the  Earth  and 
the  center  of  the  destination  planet. 

Ibe  three-dimensional  transfer  orbits  obtained  \ising  massless  departure 
and  destination  planets  are  more  realistic  than  those  obtained  using  a  simpli¬ 
fied  two-dimensional  model  of  the  solar  ^stem.  Ibe  fonaer  are  accurate 
enou^  to  give  a  good  approximation  to  the  true  values  of  the  transfer  angle, 
transfer  time,  required  velocities,  etc. 

Four  different  coeqxitatlon  procedures,  with  additional  variants  of  each, 
to  establish  the  heliocentric  transfer  orbital  elements  are  presented  and 
discussed.  The  four  different  eooQJutation  procedures  may  be  esplpyed  to 
establish  transfer  orbits  having,  respectively,  fixed  transfer  angles,  fixed 
transfer  times,  fixed  chcoacteristlc  departure  velocities,  and,  finally,  fixed 
Initial  heliocentric  velocities.  In  each  case  only  the  two-bo(^  problem  Is 
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considered  with  the  Sun  as  the  central  body.  All  of  the  methods  presented 
require  an  ephemeris  of  planets  for  the  determination  of  planeteiry  positions. 
However,  these  positions  could  be  con5)uted  accurately  enough  using  the  mean 
orbital  elements  of  the  planetary  orbits. 

The  above  requirements  are  only  a  few  of  the  more  conmon  ones  that  co\>ld 
be  considered  when  establishing  a  transfer  orbit.  None  of  the  suggested 
coBQnitational  procedures  permits  a  direct  analytical  determination  of  a 
transfer  ^ich  is  optimum  for  any  one  of  the  requirements.  TSiey  are  intended 
for  use  in  preparing  a  program  for  a  large-scale  digital  computer  idilch  would 
compute  numerous  orbits.  'Qien,  these  orbits  could  be  studied  in  order  to 
select  the  optimum  one. 


i 
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II.  GEOMETRY 

Because  the  planets  rotate  about  the  sun  in  mitually  inclined  elliptical 
orbits  and  at  different  angular  rates,  their  relative  orientations  continu¬ 
ously  change.  Uie  rate  of  change  of  orientation  of  any  two  planets  depends 
on  the  relative  size,  shape,  and  mutual  inclination  of  the  orbits.  In 
general,  for  direct  transfer  to  an  outer  planet  (larger  semimajor  axis  than 
the  departure  planet)  the  transfer  orbit  will  start  near  perihelion  for 
transfer  angles  less  than  360°.  At  the  time  of  departure  the  destination 
planet  will,  in  general,  lead  the  departiure  planet  in  its  motion  around  the 
Sun.  For  direct  transfer  to  an  inner  planet  the  destination  planet  will  lag 
and  the  transfer  orbit  will  begin  near  aphelion.  For  a  given  transfer  angle 
it  turns  out  that  the  departure  date  can  vary  over  a  relatively  short  period 
of  tiBE  without  incurring  severe  penalties.  In  other  words,  the  amount  of 
the  variation  in  the  departure  date  depends  on  the  allowable  variation  in 
the  nngnitude  and  direction  of  the  velocity  vector  of  the  vehicle  at  both 
the  departure  and  arrival  dates. 

In  the  sin^aifled  problem  which  treats  circular,  coplanar  orbits,  it  is 
possible  to  congjute  heliocentric  transfer  orbits  \islng  the  closed  analytical 
expressions  of  the  two-body  problem.  For  some  purposes,  the  orbits  established 
in  this  way  are  satisfactory  since  they  give  a  fairly  good  estimate  of  the 
departure  and  arrival  conditions.  However,  because  the  planetaiy  orbits  are 
eccentric  and  Inclined  the  orbital  data  obtained  using  the  siiiq>11fied  model 
can  differ  greatly  from  that  obtained  using  the  more  realistic  three-dimensional 
Bcdel  of  the  solar  system.  For  example,  in  the  two-dimensional  problem  the 
required  cutoff  velocity  will  decrease  as  the  transfer  angle  approaches  180°. 
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In  the  more  realistic  three-dimensional  model  the  cutoff  velocity  will  in¬ 
crease  rapidly  as  the  transfer  angle  approaches  l60°  unless  the  airrival 
planet  approaches  the  ecliptic  plane  at  the  arrival  date  (this  will  be  a 
rare  event). 
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III.  EcjUATIONS  OF  MOTIOM 


All  of  the  methods  presented  and  discussed  here  utilize  approximately 
the  same  set  of  equations  ^Ich  describe  simple  tvo-body  motion.  Only  the 
sequence  of  use  and  emphasis  of  the  various  equations  vaiy. 

The  equations  listed  below  are  in  a  general  form.  Ihe  sequence  in 
'vAilch  equations  should  be  used  and  the  syi^ols  used  are  indicated  in  each 
cethod.  3ooe  of  the  basic  equations  are  given  here  in  order  to  keep  the 
equations  and  comments  in  the  step-by-step  camputatlon  procedures  to  a  mini¬ 
mum.  %us,  the  methods  are  eSsier  to  understand  and  easier  to  cos^pare. 
Equations  based  on  two-body  motion  which  may  be  used  to  decrease  the  nxmber 
of  iterations  eire  also  given.  In  some  of  the  methods,  however,  the  use  of 
these  equations  may  require  more  coc^tation  time  than  is  saved  by  reducing 
the  number  of  iterations. 

Ihe  pertinent  equations  of  motion  and  definition  of  symbols  are  as 
follows: 

■Hie  distance  from  the  central  body  to  the  vehicle  is 


r 


a(l  -  e^) 

1  +  e  cos  (v  -  v^) 


(1) 


where 

a  >  semimajor  axis 
e  ■  eccentricl^ 

V  ■  initial  true  ancanaly  — 

o 

V  a  instantaneous  true  anomaly 

n>e  time  required  to  traverse  a  se^sent  of  an  elliptical  orbit  is. 


VO 
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accordin/];  to  Kepler's  equation 


t-t  =-  E-E  -  e  (sin  E  -  sin  E  ) 
o  n  o  0 


wliere  E  is  the  eccentric  .ononaly  and  n  is  the  mean  arqpilar  rate  in  the 
orbit.  'Ihe  cuimtities  E  arid  n  arc  obtained  frorri 


E  =  sin 


-  e‘- 


1  +  e  cos  V 


The  transfer  time  may  be  com{)utod  osin^',  Ivunbert's  Theorem.  There  are 

four  possible  elliptical  paths,  namely,  direct,  aphelion,  perilielion,  and 

(  ' ) 

indirect.  The  four  equations  vdridi  are  used  for  computing  the  transfer 
time  are 


t( direct)  =  -  q  -  -  (sin  q  -  sin 


1^11  *1 

t(aphelion)  =  - -  q  +  -  (sin  q  +  sin  q^) 


t(perihelion)  =  T  *■  q^^  -  (sin  q  +  sin  q^^^) 


?7r  i 

t( indirect)  =  - -  q  -  q^^  -  (sin  q  -  sin  q^' 


■where,  in  general 


,  /r,  +  r  + 


sin  q^/2  = 


1  /q  ^  '-2  - ' 


2  \  a 
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and 


c  = 


V 


/  2  ^  2 
r  +  r 
1  2 


2r^r2  cos 


(v  - 


V  ) 
0 


(6) 


Foi’  transfer  angles  less  than  or  equal  to  l0O°  either  Eq.  (5a)  or  (5h) 
applies.  For  anfpLes  greater  t}mn  or  equail  to  130°  Eq.  (5c)  or  (5d)  applies. 

The  transfer  ajifpLe  can  bo  exjircssed  in  terras  of  the  heliocentric  lati¬ 
tude  a,t  the  anrival  date  juid  the  difference  in  the  heliocentric  longitudes 
of  departure  and  arrival.  Tlie  equation  for  transfer  angle  is 


V  -  V 

0 


cos  \  cos  (a 


(7) 


wiiere  \  and  A  are  the  instantaneous  heliocentric  latitude  and  longitude  of 
tile  vehicle  in  tlie  heliocentric  transfer  orbit. 


IV.  OUTLIIIE  OF  R^THODS  OF  COMPUTATION 
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For  the  planninc  of  manned  or  unmanned  cpace  missions  it  will,  probably 
be  necessary  to  establish  the  transfer  orbit  subject  to  certain  constraints 
such  as  transfer  time,  depnjrture  date,  arrival  date,  guidance  and  control, 
available  thrust  enerry,  etc.  The  constraints  will  depend  in  part  on  vAiether 
tiie  mission  is  manned  or  uiimanned  and  on  the  mission  req^iireraent.  For  ex¬ 
ample,  a  manned  mission  may  require  that  the  transfer  time  be  limited  and 
specified  fairly  accurately  in  advance  of  departure.  Therefore,  it  will  be 
necessary  to  establish  a  transfer  orbit  which  will  permit  transfer  in  the 
specified  time.  Also,  if  the  mission  requires  that  the  vehicle  pass  near 
the  destination  planet  it  ma^'  be  desirable  to  either  i.unimise  the  distance 
ii-om  vehicle  to  Earth  or  to  minimise  the  velocity  of  the  vehicle  relative 
to  the  destination  planet.  Again,  this  can  be  done  if  the  transfer  orbit 
is  established  subject  to  the  proper  constraints. 

The  computation  procedures  are  presented  in  four  groups  according  to 
the  various  constraints.  V/ithin  each  group  there  are  variations  whicli  are 
indicated  by  a  quantity  called  the  piirameter.  The  groups  are  outlined  as 
follows : 

A,  Constraint  -  fixed  transfer  angle,  (v^  -  v  ,),  where  v 

13  SO.  sd 

and  v^  are  tlie  true  anomalies  of  the  vehicle  at  departure 
and  arrival 

1.  Parameter  -  t 

a 

2.  Parameter  -  t, 
d 

sd 


3.  Parameter  -  v 
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D.  Constraint  -  fiiced  tirmster  tiiae,  t 

'  c 


1.  Paratneter  -  t, 
d 


Itirai.ieter  -  t 

a 


3.  Parameter  - 


C.  Constraint  -  fixed  heliocentric  deixirture  velocity, 


1.  Ilraiietcr  -  t, 
d 

P.  Parameter  -  t 

a 

D.  Constraint  -  fixed  hyperbolic  excess  velocity, 


1.  Parameter  -  t, 
d 

V.  Parameter  -  t 

a 
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V.  PRESENTATION  AMD  DISCUSSION  OF  METHODS 


A.l  FIXED  (v  -  V  ,),  PAEIAMETER  t 

_ S _ SQ  & 

1.  Choose  t  and  obtain  A  ,  X  ,  and  r  from  an  ephemeris  or  by  compu- 

a  pa  pa 

tation.  A  and  X  are  the  heliocentric  longitude  and  latitude  of 
pa  pa 

the  destination  planet  and  r  is  the  distance  between  the  Sun  and 

pa 

destination  planet  at  the  arrival  time,  t^. 

2.  Compute  A^^  “  ^pa  "  ^^s  "  ^pa] 

3.  Assume  A^^  ”  ^sd  con^nite  r^^  from 

P 

_ e  _ 

^ed  “  1  +  e  cos  (a  ,  -  A  J 
e  ed  ep 


vhere  P  and  e  are  the  semilatus  rectiun  and  the  eccentricity  of  the 
e  e 

Earth's  ortit.  Die  longitude  of  Eajrth's  perihelion  is  A^^. 


4. 


Assume  that  r^^  «  r^^  and  compute  the  minimum  value  of  a^,  the 
semimajor  axis  of  the  transfer  ellipse,  from 


a„(min) 

8 


r  ,  + 
sd 


c 


vhere 


c 


2r  .r  cos  (v  -  v  .) 
sd  pa  s  sd' 


5.  Compute  t  from  Eq.  (5a)  if  (v  -  v  , )  <  180°  or  from  Eq.  (5c)  if 

8  8  8Q 

(v^  -  >  180°. 

6.  CoB^iute  "  ^3  (step  5)  and  obtain  A^^  and 
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7.  CcEipute  M  (step  6). 

If  M  >  0,  decrease  t^.  Use  Eq.  (5a)  or  (5c)  and  (niin)+ 

to  coii?>ute  a  new  value  of  t^.  (Note:  For  (ndn),  the  derivative 

da  /dt  -  0.  A  Ca  »  0.1  a  (min)  will  yield  a  value  for  da  /dt  which 
s'  8  8  8 

represents  approximately  the  change  cf  t^  with  a^. ) 

If  A\  <  0,  increase  t^.  Use  Eq.  (5b)  or  (5d)  and  a^  ■  a^  (min)  + 
to  coBipute  a  new  t^. 

8.  Compute  \  “  '•'a  "  7)  and  obtain  ^d  ^ed- 

9.  Compute  AA  -  Ag^  -  Ag^  (step  8). 

If  |AA|>  K  (the  size  of  K  depends  on  the  accuracy  required),  change  a^ 
(step  7)  by  an  arwunt  Aa^,  compute  t^,  and  return  to  step  8.  The  ap- 
proximte  change  in  a^  may  be  obtained  numerically  or  emalytically  from 

^  »  At  /X 

s  s' 


where 


The  miltiple  signs 
follows: 


tgOg  i  (1  -  cos  T|)  tan  n/2  -  (l  -  cos  tan  nj^/2j  (8) 
in  the  brackets  depend  on  the  type  of  transfer  path  as 


Direct  -  ♦  I%rlhellon  -  - 

i^lion  Indirect 

Ohe  eqimtion  for  X  is  derived  in  the  ^pendix. 

If  |AAj^  K,  the  transfer  orbit  is  established. 

In  A.1  the  arrival  point  is  determined  by  t^  and  the  heliocentric  longitude 
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of  the  departure  point  is  then  determined  by  using  the  coordinates  of  the 
destination  planet  at  t  and  the  transfer  angle  (v  -  v  ,). 

A  8  SQ 

Ety  assuming  that  “  ^sd  oon^jutiiig  r^^,  the  departure  point,  S^, 

of  the  transfer  orbit  is  fixed  ^en  we  assume  that  r  .  »  r 

sd  ed 

The  time  of  transfer,  t„,  is  con^wted  for  different  values  of  a  until 

8  S 

the  longitude  of  the  Eeirth  at  the  coo^juted  time  of  departure,  A^^,  is  ap¬ 
proximately  equal  to  A^^. 

Figure  1  illustrates  the  geometry  of  the  motion  for  a  typical  case  as 
the  iterations  are  performed. 


Fig.  I  —  Geometry  of  method  A.1 
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Ey  startinc  at  the  time  of  arrival,  t  ,  and  movine  backwards  in  tine, 
the  notion  of  the  vehicle  and  the  Earth  is  as  indicated  in  Fig.  1.  As  an 
exEur^e,  assiine  that  at  t  the  planet  is  at  P  and  the  Earth  is  at  E  .  During 
the  t  ,  which  is  conanited  using  the  initial  valxie  of  a„,  i.e.,  a„  (ndn),  the 

S  9 

Earth  noves  from  E^  to  E^  and  the  vehicle  noves  from  to  S^.  Ey  changing 
a  so  as  to  decrease  t  the  second  position  of  the  Earth  nay  be  E  when  the 

8  3  d 

vehicle  is  at  S  .  Additional  changes  in  a  will  cause  the  nth  position  of 

Q  3 

the  Earth,  E^^,  to  coincide  with  S^.  Uhen  this  is  accoeodlahed  the  transfer 
orbit  is  established. 


A.2  FIXED  (v  -  V  ,),  PARAMETER  t, 

_  s _ sd  ' _ ^ 

1.  Choose  t,  and  obtain  A  ,  and  r  ,,  the  heliocentric  longitude  and  distance 

d  ed  ed 

of  the  Earth  relative  to  the  Sun  at  departure. 

2.  Ckjc^iute  A  ,  the  heliocentric  longitude  of  the  destination  planet  at 

P® 

arrival,  from 


A  »  n  +  t€uj“^  (cos  1  tan  u.  ) 
pa  p  p  ^ 


where  and  i^  are,  respectively,  the  heliocentric  longitude  of  the 
ascending  node  between  the  destination  planet's  orbit  plane  cmd  the 
ecliptic  and  the  inclination  of  the  destination  planet's  orbit  plane 
to  the  ecliptic. 

The  angle  u^  «  v  is  the  angle  between  the  radius  to  the  planet 

pa  p 

at  arrival  and  the  line  of  nodes,  I.e.,  the  argument  of  latitude,  and 


It  Is  obtained  from 


u  ■  sin 


cos  (v 


v^d)  cos  P 


cot  1  cos  (o 
p  .  p 


A^d)  Bln  3 


-sin  (d  -  a  .)  cos  1  cos  3  *  cot  (v  -  v  .)  cos  (n  -  A  .) 
'  p  sd  p  '  s  sd'  '  p  sd' 
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\diere 


and 


Hie  equation  for  is  obtained  using  standard  spherical  trigonometric 
formulae  and  therefore  is  not  derived  here. 


idiere  e^  and  u)^  are  known  quantities  and  are  the  sendlatus  rectum, 
the  eccentricity,  and  the  argument  of  perihelion  of  the  destination 
planet's  orbit. 

4.  Coopute  the  minimum  value  of  a^  from 
a^  (min)  ■ 


vhere  r  .  >  r  .  and  c  is  obtained  as  in  A.l,  step  4. 
sd  ed 

5.  Compute  t^  as  in  A.l,  step 

6.  Compute  t  ■  t ,  +  t  and  obtain  A _ ,  X  ,  and  r _ . 

ads  pa'  pa  pa 

7.  Conqjute 

(v  -  V  -)  ■  cos”^  Tcos  X  cos  (a _  -  j)l 

'  s  sd'  |_  pa  'pa  sd  J 

8.  Coiniute  =.  (v^  -  v^^)(step  7)  -  (v^  -  v^^)(flxed). 

If  Av  >0.  decrease  t  :  If  Av  <0.  increase  t  # 

8  ^  s'  8  '  8 
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The  new  value  for  is  computed  as  in  A.l,  step  7. 

9.  Compute  t  =  t ,  +  t  (step  8)  and  obtain  A  .  X  ,  and  r  . 

ads 

10.  Cojmnite  (Vg  -  in  step  7. 

11.  Congnxte  Av  -  (v^  -  v^^)(8tep  10)  -  (v^  -  v^^)(fixed). 

If  |Av  |>  K,  change  a  (step  7)>  compute  t  ,  and  return  to  step  8. 
The  amount  that  a  should  be  changed,  Aa  ,  in  order  to  make  |Av  |<  K 
can  be  determined  numerically  or  analytically  by  solving  for  Aa^  as 
follows  and  returning  to  step  8. 


\4iere 


Aa^  =>  Av/X 


<i(v  -  V  -)/da„ 


(9) 


Ihe  expression  for  X,  vhlch  Is  derived  In  the  Appendix,  Is  as  follows: 


/  a  P  r 

■Y - \2  i  (1  -  cos  n)  tan  n/2  ^  (l  -  cos  tan  t]^/2 

*‘pa  L 


where  the  multiple  signs  In  the  brackets  depend  on  the  type  of  transfer 
orbit  as  follows: 


Direct  ■  ♦  Berlhelloo  ■  - 

Aphelion  ■  ♦  Indirect 

In  this  method  the  departure  point  of  the  vehicle  coincides  with  the 
Earth,  and  thtis  the  position  of  the  vehicle  at  departure  Is  determined  by 
diooslng  a  t^. 

The  first  three  steps  are  used  to  determine  the  vehicle  arrival  point 
(m  the  orbit  of  the  destination  planet,  S^.  In  general,  the  point  will 
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represent  a  past  or  future  position  of  the  destination  planet  depending  on 

the  tine  required  for  the  vehicle  to  traverse  the  specified  transfer  angle. 

The  renalnlng  stei>s  are  used  to  coiqiute  different  values  of  (v^  -  ^sd^c' 

using  values  of  t  idilcfa  depend  on  the  transfer  oxhlt^  until  the  conqiuted 

(v  -  V  natcbes  the  fixed  (v  -  v  —  In  other  vords,  until  the  po- 
'  s  sd  c  8  sd  f 

sltlen  of  the  destination  planet  coincides  with  the  position  of  the  vehicle 
at  arrival  tine. 

Figure  2  shows  the  geometiy  of  the  motion  for  a  typical  configuration 
of  the  Earth,  vehicle,  and  destination  planet.  At  the  departure  date  as- 
sune  that  the  Earth  and  vehicle  are  at  position  E^,  and  the  destination 
planet  Is  at  P..  A  point  on  the  destination  planet's  orbit,  S  ,  Is  foxuid 

d  A 

which  nakes  the  angle  between  the  radius  to  the  point,  r  ,  and  the  radius 
to  the  Earth  and  vehicle  equal  to  the  specified  (v^  -  Next,  an  a^ 
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is  chosen.  In  this  case  the  mtnlinuni  a  ,  and  the  t  is  ccniq^nited  for  transfer 
path  1.  During  the  t^  the  planet  moves  from  position  to  position  P^. 
From  Fig.  2  it  is  clear  that  the  betwen  a  radius  to  P^  and 

does  not  equal  (Vg  ~  ^sd^f*  another  t^  is  con^nxted  using  a  modi¬ 

fied  a^,  uhlch  results  in  transfer  path  number  2.  When  the  vehicle  reaches 
S  using  this  path  the  planet  msy  be  at  P  and  the  difference  between  the 

&  Ac 

specified  and  conqputed  (v^  -  v^^)  has  decreased.  Ihe  process  continues 
until  the  nth  position  of  the  planet,  P  ,  and  S  are  approximately  co- 

QA  A 

incident,  i.e.,  lAen  the  specified  and  computed  (v^  -  v^^)  are  approxi¬ 
mately  equal. 

A. 3  FIXED  (v^  -  PARAMETER 

1.  Choose  v  .. 

sd 

2.  Choose  t  and  obtain  A  ,  X  ,  and  r  . 

a  pa'  pa'  pa 

3.  CcBQute  A^^  ■  A^  -  cos“^  jcos  (v^  -  v^^)/cos  vhere  A^  is 

the  heliocentric  longitude  of  the  velocity  at  departure. 

4.  Assxnoe  that  A^^  -  A^  and  conqpute  r^  from  Bq.  (l), 

5.  Compute  a^  and  e^,  the  eccentricity  of  the  transfer  ellipse,  using 

the  equations  for  r  and  r  .,  I.e.,  Eq.  (l). 

PA  CKa 

6.  Conqxite  t^  using  Eq.  (2). 

7.  Compute  t .  ■  t  -  t  and  obtain  A  .  and  r  .. 

das  ed  ed 

8.  Ca^iute  dA  ■  A^  -  A^. 

If  |dA|  >  K,  change  t^  by  an  amount  At^  share 


Rg~208l 


At 

a 


A  (AA) 


(10) 


and  return  to  step  2. 

If  A'\  <  K,  the  transfer  orbit  is  established. 

In  the  above  procedure,  the  tine  of  arrival,  t^,  is  varied  until  the 
departure  longitudes  of  the  vehicle  and  the  Earth  are  approximately  equal. 

By  assuming  that  and  that  r^^  =  r^^,  the  position  of  the 

vehicle  at  each  departrure  date,  which  is  the  start  of  the  heliocentric 
transfer  orbit,  will  be  on  the  Earth's  orbit.  Since  (v  -  v  ,)  is  constant, 
a  change  in  the  t^  will  change  not  only  the  arrival  position  but  al.so  the 
departure  position  of  the  vehicle.  Since  the  angular  rate  of  the  Earth 
aixDund  the  Sun  differs  from  that  of  vehicle  departure  position,  a  change 
in  t  will  result  in  a  change  in  the  difference  of  the  longitudes  of  the 

cl 

Earth  and  vehicle  at  departxu'e. 

Figure  3  shows  how  the  positions  of  the  Earth  and  vehicle  at  departure 

approach  coincidence  as  the  number  of  iterations  increases.  The  angle  be- 

t-ween  the  radii  from  the  Sun  to  S  and  P  is  a  fixed  value,  namely  (v  - 

v^d)'  letters  E,  S,  and  P  denote  positions  of  the  Earth,  vehicle,  and 

the  destination  planet,  respectively. 

Initially,  the  required  position  of  the  destination  planet  at  arrival 

is  P  T .  The  vehicle  and  Earth  are  found  to  be  at  positions  3,^  and  E,^  at 
al  dl  dl 

departure. 

The  second  departure  position  of  the  vehicle,  3^,  is  computed  using 
the  second  arrival  position  of  the  destination  planet,  P  .  'Die  corresponding 
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position  of  the  Earth  at  departure  is  E^. 

OSie  tine  of  arrival  is  varied  until  the  difference  in  the  departure 
longitudes  of  the  Earth  and  vehicle  is  approximately  zero.  Assune  that 
this  condition  is  met  \iien  the  third  value  of  t  is  computed.  Then,  the 

A 

vehicle  and  Earth  will  be  at  point  E^,,  at  departure,  the  vehicle  will 
intercept  the  destination  planet  at  position  and  the  transfer  orbit 
is  established. 

The  usefulness  of  v  ,  as  a  parameter  is  not  obvious  as  in  the  case  of 
sa 

t  or  t  ;  however,  it  is  a  useful  tool  in  that  a  transfer  oihlt  can  be 
d 

established  with  a  relatively  short  con^nitation  procedure.  Also,  the  ele¬ 
vation  aivgle  of  the  vehicle’s  velocity  vector  at  departure,  7^^,  changes 
directly  with  v^^.  In  general,  as  v^^  increases,  7^^  will  increase,  and 
as  a  consequence,  required  hyperbolic  excess  or  cutoff  velocity  will  in- 
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B.l  FIXED  t  ,  PARAMETER  t, 

_ 3*  Q. 

1.  Choose  and  obtain  and 

2.  Compute  *  ^d  ^  obtain  A^^,  X^,  and 

3.  Compute  (Vg  -  from 

^''s  '  ""sd^  “  [°°®  ’^pa  ^'^pa  •  ^sd^] 

^ere 


A  j  -  A  , 

sd  ed 


4.  f-nmpit/*  the  wi-i  ni  imm  valvie  Of  a  from 


a^  (min) 


r  j  +  r _  +  c 

sd  pa 


idiere  r  ^  =  r  and  c  is  obtained  from  Eq.  (6). 
sd  ed 

3.  Compute  t^  as  in  A.l,  step 

6.  Cosqinite  Ats  -  t^  (step  5)  -  (fixed). 

If  At  /  0,  conmute  a  new  t„  as  in  A.1,  step  7« 

8  S 

7.  Conqnite  At^  -  t^  (step  6)  -  t^  (fixed). 

If  |Atg|  >  K,  change  a^  (step  6)  by  an  amount  Aa^,  camnite  t^  and  At^. 

Repeat  step  7  until  the  orbit  is  established,  i.e.,  until  <  K. 

The  amount  of  change  for  a  nay  be  coimnited  as  in  A.l,  step  9« 

With  t^  a  given  constant,  a  choice  of  t^  determines  the  required  t^. 

The  heliocentric  coordinates  of  the  Earth  at  t^  and  those  of  the  destination 
planet  at  t^  represent  two  fixed  points  separated  by  a  constant  (v^  -  v^^) 
which  is  determined  by  the  coordinates. 


RH-2681 

21 


In  order  for  the  vehicle  to  Intercept  the  destination  planet  at  t^, 
the  actual  t  ,  \diich  depends  on  the  transfer  orbit  used,  mist  eq^  the 
given  tg.  Ihis  is  accomplished  by  changing  a^,  uhlch  causes  a  change  in 
the  actual  t^. 

In  Fig.  U  the  point  E^,  represents  the  position  of  the  Earth  and 
vehicle  at  t,.  T5ie  point  S  represents  the  position  of  the  destination 
planet  at  the  required  t^  which  is  fixed  by  the  given  t^  and  chosen  t^. 

AssuiiB  that  the  semlmajor  axis  of  transfer  orbit  number  1  is  a(min) 
and  that  the  corresponding  t  is  larger  than  the  given  t  .  'Dien,  it  is 
that  the  vehicle  will  arrive  at  position  S  too  late  and  the  desti- 

cL 

nation  planet  will  have  moved  to  position  P^.  If  the  correct  equation 

for  t  is  used  with  an  increased  a  ,  the  actual  t  will  decrease  and  the 
s  s  * 
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vehicle  ney  arrive  on  transfer  orbit  number  2  either  too  early  or  too  late. 

According  to  Fig.  4,  arrival  la  too  late,  l.e.,  t^  Is  too  large,  and  the 

destination  planet  Is  at  P  >*en  the  vehicle  Is  at  S  .  Further  changes 

cud 

In  a  will  result  In  a  t  ^Ich  will  permit  the  vehicle  to  travel  on  transfer 

S  8 

orbit  nuinber  n  and  Intercept  the  destination  planet  at  S^. 

B.2  FIXED  tg,  PARAMETER  t^ 

1.  Choose  the  t^  and  obtain  and  r^. 

2.  Ccn?)ute  ^  ~  obtain  and  r^^. 

These  steps  are  the  same  as  In  B.l. 

Since  this  method  differs  only  slightly  from  B.l  ■vdiere  t^  was  the 

paraoeter,  the  discussion  of  the  method  Is  similar  to  that  of  B.l,  and 

Fig.  4  serves  to  Illustrate  the  geometry. 

A  <dioiee  of  the  t  fixes  the  rendezvous  position  of  the  vehicle  and 
a 

destination  planet  at  point  S^.  Since  the  t^  Is  fixed,  a  choice  of  t^ 
determines  uniquely  the  required  t^.  Thus,  the  required  departure  position 
of  the  vehicle,  which  must  coincide  with  the  Earth's  position  at  departure. 

Is  determined. 

As  in  B.l,  a  Is  varied  until  the  cooiinxted  t„  matches  the  given  t  . 

'  B  8  8 

Uhen  this  occurs,  the  Earth  and  vehicle  are  at  the  point  E^,  at  the 
required  departure  date  and  the  transfer  orbit  Is  established. 

B.3  FIXED  t^,  PARAMETER  v^ 

1.  Choose  v^. 


^  -  ^d  ^  % 


2.  Choose  t^  and  conqute 
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3.  Obtain  r^^.and  for  and  t^,  respectively. 

4.  Coapute  v  from 


sd 


+  cos”^  cos  cos  (a 


pa 


pa 


-''.a)] 


idiere 


A  ,  ■  A  , 
sd  ed 


5.  CoB^Jute  e„  and  a  using  Eq.  (l). 

S  8 

6.  Coiqnite  the  minimum  value  of  a  from 


S(min)  -  i  +  c) 


vbere  r^^  =  r^^  and  c  is  obtained  using  Eq.  (6). 

7.  CoBQaite  t^  for  a^Cmin)  \ising  Eq.  (^)  or  (^c)  depending  on  the  size 

^^s  - 

8.  CoB^te  Atg  =  tg  (step  7)  -  (fixed). 

If  At  /  0  xise  either  Eq.  (5a)  or  (5b),  depending  on  the  "^pe  of 
s 

transfer  orbit  desired,  vith  a  (step  5)  to  compute  t  if(v  -v.) 

8  8  8  8Q 

^  180°.  Similarly,  if  (v  -  v  ,)  >  l80°,  use  either  Eq.  (5c)  or  (5d). 


If  At  >  0,  decrease  t  (step  7)  by  using  a„  (step  5)  with  Eq.  (5a) 
8  8  8 

or  (5c)  depending  on  the  size  of  (v^  -  Vg^)» 

If  At  <  0,  increase  t  (step  7)  by  \islng  a.  (step  5)  with  Eq.  (5b) 
8  8  8 

or  (5d)  depending  on  the  size  of  (v^  -  ^sd^* 

Caequte  At^  ■  t^  (step  8)  -  t^  (fixed). 

If  |At^|  >  K,  diange  t^  by  an  asiount  At^  and  return  to  step  2. 


9. 


ra'i-208l 


At^  =  At  /X 
d  s' 


vhere 


.  /a  P  , 

X  I 

r 

pa 


,“s<^  *  =s)  -  ■'sd 


sin  V 


(l  -  e  )(r  -  r  ) 

s  sd  pa 


^  2  "n 

^  -  (l  -  cos  T])  tan  ^ 


-  (1  -  cos  T]^)  tan  -— 


sd  .  /  \ 

+  x; —  sin  (v  -  V  ,) 

2c  ^  s  sd' 


+  T1  +  '1 

-  tan  ^  -  tan  — 


(11) 


The  multiple  sipiits  in  the  brackets  depend  on  the  type  of  transfer  path 
as  follows: 


Direct  =  -,  +  and  +,  +  Perihelion  -  and  +,  - 

Aphelion  =  +,  +  and  +  Indirect  -  -  and  - 

The  equation  for  X  is  deri’/ed  in  the  Appendix. 

The  discussion  concerning  the  usefulness  of  as  a  parameter  is 
given  in  method  A. 3. 

The  geometry  of  this  method  is  illustrated  by  Fig.  5- 
A  selection  of  t^  gives  a  position  for  the  Earth  at  departure  E^  and 
because  t,  is  fixed,  t  is  determined,  and  consequently  the  cooixiinates  of 

S  3. 

the  planet  position,  P^,  at  the  arrival  time  are  known. 

The  time  required  for  the  vehicle  to  travel  from  E^  to  P^  is  computed 
using  a  prescribed  v^^.  If  the  computed  t^  differs  from  the  fixed  t^,  the 
vehicle  will  not  intercept  the  plam^t  at  P^.  Since  the  t^  will,  in  this 
case,  vary  directly  with  (v^  -  a  cliange  in  (v^  --  v^^)  will  cause  a 

change  in  t^. 
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A  second  value  of  t^,  and  consequently  vd.ll  fix  the  Earth  at  point 
and  the  planet  at  ^ince  the  angular  rates  of  the  planets  around 

the  Sun  are  different,  the  angle  between  the  radii  to  points  E^  and 
will  not  equal  the  ang].e  between  the  radii  to  points  E^  and  P^.  Conse¬ 
quently,  tlie  t^  vd.ll  be  different, 

J]y  varying  t^  the  correct  i-elative  orientation  of  die  deivarture  .and 
destination  planets  can  be  found  vihich  vd.U.  give  the  correct  t^. 


Fig.  5 — Geometry  of  method  B.3 
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C.l  FIXED  PARAMETER 


1.  Choose  t,  and  obtain  A  ,  and  r  .. 

d  ed  ed 


2.  Compute  a  from 
s 


kA  r  j 

8  sd 


■  ^sd'^sd 


^d  “  ’^ed* 


3>  Choose  t  and  obtain  A  ,  X  ,  and  r  . 

a  pa'  pa'  pa 

1*.  Coomute  t  =  t  -  t,. 

sad 

5.  Compute  (v  -  v  ,)  from 

s  Sd 


(v  -  V  , )  ■  cos"^  fcos  X  cos  (a  -  a  .)! 
3  sd'  [_  pa  'pa  sd'J 


6.  Compute  the  minimum  value  of  a  from 

s 


a^(mln)  -  1/4  (^3^  +  ♦  c) 

vhere  c  Is  obtained  from  Eq.  (6). 

7.  Congiute  Aa  ■  a  (min)  -  a  (step  2). 

8  S  S 

If  Aa  >  0,  change  t  by  At  idiere 

8  Ci  & 

Ata  -  AaaA 


^here 


sin  (a  -  a  .) 
'  pa  sd' 


X 


(12) 
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and  return  to  step  3.  The  equation  for  X  is  dervied  in  the  Appendix. 

If  ^  0,  proceed  to  step  8. 

8.  Compute  t  for  a  (min)  using  Eq.  (5a)  or  (5c)  depending  on  the  size 

S  8 

-  '^sd)* 

9.  Compute  At_  t  (step  8)  -  t  (step  4). 

3  8  3 

If  At  =  0,  either  Eq.  (5a)  or  (5l>),  depending  on  the  type  of  transfer 
8 

orbit  desired,  can  be  used  with  a^  (step  2)  to  conipute  t^  if  (v^  -  v^^) 
^180°.  Similarly,  if  (v^  -  v^^)  >  180°,  either  Eq.  (5c)  or  (5d)  would 


be  used. 


If  At  >  0,  decrease  t„  (step  8)  by  using  a  (step  2)  with  Eq.  (5a) 
or  (5c)  depending  on  the  size  of  (v^  -  v^^). 

If  At  <0,  increase  t  (step  8)  by  using  a  (step  2)  with  Eq.  (5b) 

or  (5d)  depending  on  the  size  of  (v^  -  v^^). 

10.  Compute  At  ■  t  (step  9)  ■  t  (step  4). 

8  8  8 

If  lAt  I  >  K,  change  t  by  At  and  return  to  step  3. 

At  ■  At  /X 
a  s' 


\diere 


gc  f  sin  (Ap^  -  Ag^)(-  tan  n/2  i  tan  ^^/2)  -  1 


(13) 


pa 


'nie  multiple  signs  in  the  parentheses  depend  the  type  of  transfer 
path  &a  follows: 

Direct  ■  +,  +  Perihelion  ■  - 


Aphelion 


Indirect 


The  equation  for  X  is  derived  in  the  Appendix. 
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If  |Atg|  <  K  the  transfer  orbit  is  established. 

As  in  A.l,  step  9i  the  equation  for  X  will  depend  on  the  size  of 
(v^  -  and  the  type  of  transfer  path. 

This  computation  procedure  may  be  used  to  obtain  transfer  orbits  for 
given  heliocentric  departure  velocities  for  a  specified  departure  date. 

Ety  using  t^  as  a  parameter  and  assuming  that  the  IJarth  and  vehicle  are 
coincident  at  departiue,  a  search  for  the  correct  trajectory  is  made  by 
vaiylng  t^  until  the  vehicle  intercepts  the  destination  planet.  For 
and  t^  some  arrival  times  cannot  be  realized,  i.e.,  the  a^  vdiich  is  de¬ 
termined  by  and  t^  will  be  less  than  the  required  rndnimum  which  takes 
into  account  the  distance  to  the  destination  planet  at  arrival  and  the 
transfer  angle,  'nius,  a  change  in  t  will  be  required. 

Figure  6  shows  the  geonetry  of  motion  for  a  typical  case  which  utilizes 
this  cosputation  procedure. 

The  departure  position  of  the  Earth  and  vehicle,  E^,  S^,  is  determined 
by  a  (dioice  of  t^.  The  time  of  arrival  is  arbitrarily  chosen  and  the  helio¬ 
centric  coordinates  of  the  destination  planet  are  obtained.  Assume  that 
the  planet  is  at  point  S  at  the  time  of  arrival.  Ihe  coordinates  of  points 
and  are  used  to  cospute  the  minimum  value  of  a^.  This  pezmlts  one 
to  determine  if  the  given  a  Is  large  enough  to  yield  a  transfer  orbit. 

If  the  given  a^  is  less  than  the  minimum  a^  (see  step  T)  &  new  value  of  t^ 

Is  determined.  Otherwise  at  is  conputed,  using  the  minimum  a  ,  which 
corresponds  to  transfer  orbit  nunter  one.  If  the  is  too  large  the  planet 
may  be  at  whan  the  vehicle  arrives  at  S^.  The  relative  position  of 
and  permits  one  to  select  the  correct  equation  (see  Eq.  9)  for  t^  when 
transfer  orbit  number  two  Is  determined  using  the  given  a^. 
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In  general  the  first  valvie  of  t  will  not  be  correct  and  several  values 

& 

of  t  will  be  used  before  the  correct  transfer  orbit  Is  established, 
a 


Fig.  6 — Geometry  of  method  C.1 
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C.2  FIXED  PARAMETER 

1.  Choose  t_  and  obtain  A  ,  \  ,  and  r  . 

a  pa  pa  pa 

2.  Choose  and  obtain  and 

3.  Conpute  tg  »  t^  -  t^. 

4.  Conpute  a^  from 


4  r, 


sd 


-  >^sd^sd 


^d  “  'ed- 


5.  Conpute  (Vg  -  from 


(v  -  V  j)  »  cos 
'  s  sd' 


-1 


cos  X  cos  (a__  -  A  ,) 
pa  P®  3®  I 


6.  Conpute  the  odniimim  value  of  a^  from 


a^(inin)  -  J  +  V  ^ 

^ere  c  is  obtained  using  Eq.  (6). 

7.  Compute  =»  a_(inin)  -  a  (step  4). 

6  8  8 

If  oa,^  >  0,  change  t^  by  an  amount  At^  and  return  to  step  2. 

At.  ■  Aa_/X 

d  8 


vhere 


X  - 


4c  r 


sin  (a 


sd 


-  V 


(14) 


The  derivation  of  the  equation  for  X  is  similar  to  the  derivation  of 
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the  equation  for  X  in  C.l,  step  7  and  therefore  is  not  included  in 
the  Appendix. 

If  Aa^  <  0,  proceed  to  step  8. 

8.  (loinpute  t  for  a  (min)  using  Eq.  (5a)  or  (5c)  depending  on  the  si?e 

s  s 

of  (v  -  V  J. 

S  sd 

9.  Compute  At  =  t  (step  8)  -  t  (step  3)- 

s  s  s 

If  At^  =  0  either  Eq.  (5a)  or  (5h);  depending  on  the  type  of  transfer 
orbit  (desired,  can  be  used  ’.hth  a  (step  4)  to  compute  t  if  (v  -  v  ,) 

S  3  S  SGL 

<180°.  Ciriilarly,  if  (v  -  )  >180*^,  either  Eq.  (5c)  or  (5d)  can 

be  used. 

If  At  >  0,  decrease  t  (step  8)  by  using  a  (step  4)  with  Eq.  (5a) 

s  s  s 

or  (5c)  depending'  on  the  size  of  (v^  -  v^^). 

If  At  <  0,  increase  t  (step  8)  by  using  a  (step  4)  with  Eq.  (5b)  or 

s  s  s 

(5d)  depending  on  the  size  of  (v^  - 

10.  Compute  At  =  t  (step  9)  -  t  (step  3)» 

S  6  S 

If  |At^|  >  K,  change  t^  by  At^  and  return  to  step  2. 

At,  =  At  /X 
d  s' 

where 


a  P  r 
set 


2r  ,c 
sd 


sin  [h  -  A  .)  -  tan 


The  signs  of  the  bracketed  terms  ore  determined  as  in  C.l,  step  10. 
The  derivation  of  the  equation  for  X  is  similar  to  tliat  for  At^  in 

C.l,  step  10,  and  therefore  is  not  given. 

If  |At^|  <  K,  the  transfer  orbit  is  established. 
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Ghis  method  of  establishing  a  transfer  orbit  is  similar  to  C.l.  13ie 
major  difference  is  that  the  time  of  arrival  is  fixed  and  the  departvire 
date  is  used  as  a  variable.  Becaiise  t^  is  varied  it  is  necessary  to  com¬ 
pute  a  new  value  of  a^  for  each  iteration,  \Aiile  in  C.l  a  choice  of  t^ 
fixed  the  value  of  a^  for  all  the  cos^tations. 

Die  time  of  arrival  is  chosen  and  we  assume  that  the  vehicle  is  coin¬ 
cident  with  the  destination  planets  at  arrival.  Next  t^  is  varied  until 
a  transfer  orbit  is  found  \diich  will  make  the  departure  coordinates  of  the 
Earth  and  vehicle  equal. 

Similar  to  C.l,  some  departure  dates  cannot  be  used  for  the  given 
and  selected  t  .  Ihat  is,  the  required  V  or  a  may  exceed  the  given 
values.  In  this  event  it  is  necessary  to  select  another  departure  date  so 
as  to  reduce  the  required  a^. 

Figure  7  shows  the  geooetry  and  relative  positions  of  the  planets  and 

vehicle  using  this  coiqmtation  procedure. 

Ihe  geometry  of  the  notion  is  easy  to  understand  if  we  start  at  t  and 

move  the  planets  and  vehicle  backwards  in  time.  At  t  the  destination 

& 

planet  and  vehicle  are  at  point  S  ,  P  on  the  destination  planet's  orbit 

Cl  Cl 

and  the  Earth  is  at  point  E  .  Assume  that  transfer  orbit  1  corresponds 
to  the  mini  nan  value  of  a  .  During  the  corresponding  time  t  the  vehicle 

B  S 

moves  to  point  while  the  Earth  may  move  to  point  E^.  In  that  case  the 

t  is  too  large  and  consequently  the  caqnited  t.  is  too  tnnali  (early), 
s  d 

Ihe  relative  positions  of  the  Earth  and  vehicle  at  departure,  l.e.,  E^ 
and  indicate  whldi  part  of  Eq.  (3)  should  be  used  with  the  given  a^  to 
eoBvute  a  new  value  of  t  .  If  orbit  2  correspcatds  to  the  given  a  ,  then 

B  B 

the  Earth  ney  be  at  E.^  vhen  the  vehicle  is  at  S,.  Since  B.^  and  S.  are 
02  d  02  o 


D.l  FIXED  V  ,  PARAf-ETER  t, 
a/ _ d 

1.  Choose  t,  and  obtain  A  ,  and  r 

d  ed  ed 
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2.  Choose  t  . 

s 

3.  CoB^Jute  \  *3  obtain  and 

4.  Compute  (v^  -  =  cos"^  jcos  cos  (a^^  -  A^^) 

^ere 


A  ,  =  A  , 
sd  ed 


5.  Choose  and  coignite  from 


'^s  “  ^^s  *  5)  + 

6.  Compute  e  and  a  using  Eq.  (l). 

8  S 

7.  Coiqpute  the  minimum  a^  from 

a^(min)  -  J  (r^^  * 

idiere  r  .  =  r  ,  and  c  Is  obtained  from  Eq.  (6). 

Sd  eel 

8.  Compute  t  using  a  (ndn)  and  Eq.  (3s)  or  (3c)  depending  on  the  size 

S  3 

^^8  -  ^Sd^- 

9.  Compute  At  at  (step  8)  -  t  (step  2)  and  compute  t  as  In  method 

S  S  8  S 

C.2,  step  9- 

10.  CoB?>ute  At  -  t  (step  9)  -  t  (step  2). 

8  8  8 

If  |Atg|  >  K,  change  (step  3)  ^  ^^sd  return  to  step  5* 


Av  ,  a  At  /X 
sd  s' 


^ere 
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X  = 


s  1  -  e 


,(l  +  ef)  -  r 


sd 


r  .  -  r 
sd  pa 


(r  sin  V  -  r  ,  sin  v  ,) 
pa  s  sd  sd 


r  ,  sin  V  , 
sd  sd 


+  n  t  -  (l  -  cos  n)  tan  T1/2 
2  s  s  ' 


-  (1  -  cos  tan  1)^/2 


(16) 


The  niultiple  signs  in  the  brackets  depend  on  the  type  of  transfer  path 
as  follows: 

Direct  =  -,  +  Perihelion  =  - 

Aphelion  =  +,  +  Indirect  =  +,  - 

Tlie  equation  for  X  is  deri'/ed  in  the  Appendix. 

If  [At^l  <  K,  proceed  to  step  11. 

11.  Compute  the  Earth's  velocity  vector  at  departure  from 


V 


ed 


/2  1  X 

^^s  (7-  -  7") 

ed  e 


where  a^  is  the  seraimajor  axis  of  the  Earth's  orbit. 
12,  Compute  the  elevation  angle  of  from 


7  ,  =  cos 
ed 


-1 


2  2 \ 

,  ,  a  (1  -  e  ) 

+  /  e  e 

r  ,  (2a  -  r  ,) 

ed  e  ed' 


where  the  plus  sign  is  used  if  r^^  >  0. 


13.  Compute  the  TOhicle's  velocity  vector  at  departure  fixim 


sd 


^2  1  ^ 
sd  s 
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14.  Compute  the  elevation  angle  of  from 


lAere  the  plus  sign  is  used  if  0  <  <  n. 

15.  Compute  i  ,  the  inclination  angle  of  the  transfer  orbit  plane,  from 

s 

i  =  tan  ^  tan  X  /sin  (a  -  A  J 
s  pa'  “pa  sd 

16.  Compute  0:,  the  ajigle  bet'v.’een  and  from 


Note:  no  attempt  is  made  here  to  express  X  in  terms  of  the  orbital 
parameters  since  the  equation  is  tedious  even  with  several  simplifying 
assumptions.  A  purely  iterative  procedure  is  suggested. 
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If  IaV  I  <  K,  the  orbit  is  established. 

I  co’ 

In  order  for  a  transfer  orbit  to  be  possible,  the  specified  must  be 
equal  to  or  greater  than  the  mlnirainn  value  of  V^.  For  transfer  to  inner 
or  outer  planets,  a  >  ^  (r  +  r  )  -where  the  r's  ore  the  perihelion 
dls-tances  of  -tlie  orbits  of  Barth  and  the  destination  planet. 

For  transfer  to  outer  planets,  must  satisfy  the  equation 


V 

00 


'V 


r  /r 
pp'  ep 


PP 


r  ) 

ep' 


For  trojisfer  to  inner  planets 


a 

e 


) 


V  >  /pji  r  /r  (r  -i-  r  )  -  /ii- - ■^) 

oo  “  '\y  pa'  ea  ea  pa  -yy  r^^  a^ 

■vdiere  r  and  r  are  -Uie  aphelion  distances  of  the  planet  and  Earth. 

Ihe  geometry  of  the  motion  for  a  special  transfer  is  sho-wn  in  Fig,  8. 

The  first  -ten  steps  of  this  method  are  used  to  establish  an  orbit  be- 
t-ween  a  known  departure  point  and  a  known  arri-val  point,  which  i-fill  pemit 
the  -vehicle  to  intercept  -the  destination  planet.  Hie  remaining  steps  of 
■the  me-thods  are  used  to  determine  the  relati-'re  magnitudes  of  -the  computed 

-which  is  based  on  the  es-tablished  transfer  orbit,  and  the  specified  V^. 

Hie  choice  of  the  parame-ter  t^  determines  tlie  position  on  the  Earth's 

orbit  E^,  3^.  This  point  remains  fixed  -vdiile  v^^  and  then,  if  necessary, 

t  is  varied  in  order  to  es-tabllsh  a  trauisfer  orbit, 
s 

A  selection  of  t^  and  t^  determines  the  positions  E^,  S^,  and  S^,  -which 
is  a  position  on  -the  destination  planet's  orbit.  Assume  -that  orbit  1  is 
computed  using  a  (rain)  and  that  t  is  too  large.  In  -this  case  the  planet 
may  be  at  P  ^  -when  the  vehicle  is  at  3  .  Orbit  2  is  computed  using  -the 
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a  •vAiich  is  based  on  the  coordinates  of  S,  and  S  ^  and  the  v  The  t 
s  d  ai  sd  B 

may  be  su'Ch  that  the  planet  is  at  vehicle  is  at  S^.  At  this 

point  in  the  computations  the  is  changed  and  new  orbits  are  computed 
until  the  orbit  n  causes  P  and  S  ,  to  coincide,  i.e.,  the  vehicle  inter- 
cepts  the  planet. 

For  orbit  n,  which  connects  S,  and  S  the  required  V  is  determined 
and  compared  to  the  specified  V^.  If  the  V^'s  disagree,  then  the  arrival 
point  on  the  destination  planet  orbit  is  changed  by  taking  a  different 
value  of  t^.  As  before,  orbits  are  computed  until  the  vehicle  intercepts 
the  planet  at  S  .  The  V  's  are  again  compared,  and  if  they  are  different 
a  new  value  of  t^  is  determined  and  the  process  repeats  until  an  orbit  is 
foiuid  wdiich  will  cause  the  reqxiired  to  match  the  specified  V^.  If  no 
such  orbit  can  be  found  then  the  parameter  t^  must  be  changed. 

D.2  FIXED  V  ,  PARAMETER  t 
_ _ a 

1.  Choose  t  and  obtain  A  ,  X  ,  and  r  . 

a  pa'  pa'  pa 

2.  Choose  t  . 

8 

3.  Compute  t ,  =  t  -  t  and  obtain  A  ,  and  r 

das  ed  ed 

1<--18.  These  steps  are  the  same  as  in  D.l. 

The  geometry  of  motion  is  different  from  that  of  D.l  and  is  shown  in 
Fig.  9. 

In  this  method  the  tijos  of  arrival  is  selected,  consequently  the  position 
of  the  destination  planet  at  the  arrival  time.  This  method  is  similar  to 
D.l  in  that  v^^  is  used  as  the  variable  in  computing  the  diffei'ent  orbits 
which  connect  the  fixed  arrival  point  and  a  particular  departiLre  point  on 
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the  Eaiih's  orbit;  it  differs  from  D,1  in  that  a  change  in  t  causes  the 

s 

departure  point  to  move  along  the  Earth ' s  orbit  while  the  arrival  point 
remains  fixed. 


Fig.  9  —  Geometry  of  method  D.2 
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VI.  DISCUSSION 


Die  methods  presented  and  discussed  here  are  only  a  few  of  many  which 
could  be  employed  to  establish  interplanetary  transfer  orbits.  These 
methods  are  not  necessarily  better  than  others,  but  they  represent  some  of 
the  more  obvious  approaches  to  the  problem  and  at  the  some  time  involve 
some  of  the  parameters  of  greatest  interest. 

The  detailed  computation  procedure  given  with  each  method  will  be  of 
value  in  preparing  a  computer  program.  The  computation  procedure  plus  the 
discussion  and  figures  will  give  the  reader  a  "feel"  for  the  relativne  ge¬ 
ometry  of  motion  and  spatial  orientation  of  the  vehicle  and  planets  as  a 
function  of  time  for  each  of  the  various  parameters. 

Several  equations  ore  included  in  the  computation  procedures,  namely 
Eqs.  (8)  -  (l6),  for  the  pxirpose  of  reducing  the  number  of  iterations  re¬ 
quired.  These  equations  are  obtained  by  taking  the  derivatives  of  various 
equations  of  two-body  motion.  Since  the  evaluation  of  these  equations  is 
straightforward,  though  perhaps  tedious,  their  use  may  in  some  cases  reduce 
the  computation  time  that  would  be  required  using  a  purely  iterative  pro¬ 
cedure.  The  advantage  of  using  these  equations  rather  than  a  simple  trial- 
and-error  approach  depends  soraevdiat  on  the  type  of  computer  used. 
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Appendix 

DERIVATION  OF  E.^UATIONS 

DERIVATION  OF  (8)  OF  METHOD  A.l 

Acconiing  to  Eq.  {5),  t  is  a  function  of  n  ^  t],  nnd  t]  .  Thus 

S  S  J_ 


where 


dt„ 

da 


17^  ~  „3/2 


and  wtiere 


>1=2  sin 


-1 


=  2  sin 


-1 


,  ,  r  ,  +  r  +  c 
X  j  sd  pa 

2  . 

B 


T  /  r  ,  +  r  -  c 
1/  sd  pa 

2 


1/2 


1/2 


c  = 


2 

r  ,  +  r  -  2r  ,r  cos  (v  ■■  v  , ) 
sd  pa  sd  pa  s  sd 


1/2 


(17) 


(18) 


(19) 


(20) 


(21) 


For  departures  from  Earth  and  heliocentric  transfer  orbits,  r^^  =  r^^  and 

P  is  the  heliocentric  gravitational  constant.  For  con\nsnience  the  units 
8 


=  1. 


may  be  chosen  so  that  k 
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The  total  derivatives  are  evaluated  as  follows 


d  — - 


^'1  ^  ^  ^  ^  +  ^ii_ 

da  8r  ,  da  or  da  oc  da  oa 

s  sd  s  pa  s  s  s 


Since  t  is  the  parameter,  r  is  constant.  V/e  assume  that  r  ,  is  a 
a  pa  sd 

constant. 

Also,  because  (v  -  v  .)  is  fixed,  c  is  constant.  Tlius 
'  3  sd  ’ 


=  -  —  tan  r]/2 

da  oa  a  ' 


In  a  similar  manner  we  find  that 


tan  Ti^/2 


The  partial  derivatives  of  t^  depend  on  the  type  of  transfer  orbit  and 
are  obtained  using  Eq.  (5).  The  expressions  for  the  partiols  are 


8t  t 

S  H 

S  £ 


for  all  types  of  transfer  orbits. 


3  +1 


(1  -  COB  q) 


(26) 
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'vAiere  the  plus  sign  is  used  for  direct  and  perihelion  transfer  aiid  the 
minus  sign  is  used  for  aphelion  and  Indirect  transfer.' 

S  _  i  ^  (1  .  cos  Til) 

X  8 

Tdiere  the  plus  sign  is  used  for  perihelion  and  indirect  transfer  and  the 
minus  sign  is  used  for  direct  and  aphelion  transfer. 

Substitution  of  the  derivatives  into  Eq.  (17)  gives 

•|  t^n^  -  (1  -  cos  T|)  tan  T1/2  -  (l  -  cos  tan  ti^/2 

(26) 

vfaere  the  multiple  signs  in  the  brackets  depend  on  the  type  of  tremsfer  as 
follows: 

Direct  ■  -,  +  Bsplhellon  ■  - 

Aphelion  ■  ♦  Indirect  •  - 


DBRVIAIION  OF  Evi.  (9)  OF  MBTOOD  A.2 

d(^  -  ^d>  -  -  ^d^  -  ^d^ 

”  ^s  "  ^s 


(29) 


(v^  -  v^)  is  related  to  the  coordinates  of  the  departure  and  destination 
planets  as  fallows: 


V 

8 


8d 


■  cos”^  ^cos  cos  (a 


pa 


pa 


(30) 


Ittt-2681 

1*6 


'^sd  ’  '^ed* 


-  V)  -  "sd)  ^  ^  -  ’'.d)  ^  (3,) 

—s; - "p. 

Since  t^  is  the  parajneter  in  this  method,  le  constant  and 


dA 


sd 


da 


Now,  the  total  derivatives  of  Eq.  (3I)  are 
SX  dt 


dX 

— 2S  ,  ^ 

da  ot  da 
s  a  s 


■  55“ 


da 


(32) 


Since  t^  is  fixed 


JBS 

da 

s 


Sx_  dt 
pa  _ B 

dt  da_ 
a  s 


Similaply 


_SS 

da 

s  a 


^  IT 


According  to  Eq.  (5),  ^3  1“  »  function  of  n^,  n,  and  idiich  are 
given  hy  Eqs.  (I8)  -  (20). 

A  (dioice  of  t,  fixes  r  ..  We  assume  that  r  is  constant;  since 
a  SCI 

(v  -  V  )  is  fixed,  c  is  constant  and  dt  /da  is  given  ty  Eq.  (26), 

\  g  '  SO 
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Is 


da 


V  ^ 


I  -  (1  -  cos  n)  tan  11/2  -  (l  -  cos  tan  T)j^/2j  (28) 


Ihe  signs  of  the  terms  are  the  sane  as  in  the  derivation  of  Eq.  (8). 
Now,  from  the  geonetiy 


X  =  tan~^  I  t€in  1  cos  (A _ 

pa  L  p  '  pa 


■v>] 


(33) 


and 


5x  Sa 


3^  -  3^  tan  ip  cos  Xp^  cos  (Ap^  - 


(34) 


The  partial  derivative 


m  A 


pa 


pa 


(35) 


and  Eq.  (34)  becomes 


St^  "  V  ”  ~~"2  ^  V 


ISie  partial  derivatives  of  (v^  -  v^^}  are  foimd  \ising  Eq.  (30).  ’Qiey 


^^^s  '  ^sd^  Bin  Xp^  cos  (A^  -  A^) 

3x—  -  ein  (v^  -  vj^) 


(37) 
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and 


8(v  -  V  J 

'  s  sd' 


IT 


p<i 


cos  k  sin  (a  -  a  ,) 
_ pa  ^  pa  sd' 

sin  (v  -  V  ,) 

'  a  sd 


(38) 


Substitution  of  Eqs^  (3^0,  (35)>  (37);  and  (38)  into  Eq,  (3I),  after 

usin^^  flTml  1  -angle  approximations,  i.e.,  sin  k_^^  =  cos  =  1,  and 

tan  i  =  i  ,  rives 
P  P 


d(v  -  V  ) 
s  sd 

da 


r  ^  sin  (v  -  v  , ) 
pa  s  sd 


k  i  cos  (a  -  a  ,) 
pa  p  pa  sd 


dt 

c 

da 


Since  and  i^  are  small  for  the  planets  of  interest,  their  product 

is  neglected.  In  accord  with  this  and  the  assumptions  above 

sin  (v  -  V  ,)  I  sin  (A  -  A  ,) 
s  sd  '  pa  sd' 


and  we  get 


d(v  -  V  J  ^p  P  dt 

'  f!  BCi  V  s  p  8 


3 _ sa 

da 


2  da 
r  s 

pa 


Substituting  for  dt^/da  gives 
s  s 


d(v  -  V  ,) 
s  sd 

da 


a  P 


pa 


+  ^  t  n  -  {: 
2  s  3 


-  cos  n)  tan  q/2  -  (1  -  cos  q,  )  tan  q-/2 


(39) 
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where  the  nmltiple  signs  in  the  brackets  depend  on  the  type  cf  transfer 
path  as  follows : 

Direct  =  ^  i'erihelion  -■  - 

Aphelion  =  -*,  +  Indirect  =  +,  - 

DERIVATION  OF  Eg.  (lO)  OF  t^ETHOD  A.  1 


d(AA  )  _  ^'ed^ 


dA 


sd 


dA 


ed 


dt 


dt 


dt 


dt 


(40) 


Since  (v^  -  is  fixed  and  is  specified^  is  determined.  If,  in 


sd 


addition,  we  assume  that  r  ,  and  r  are  constant,  t  will  also  be  constant, 
’  sd  pa  ^  s  ’ 

and  from  t  =  t,  +  t 
ads 


dt  =  dt^ 
a  d 


Eq.  (40)  becomes 


dAA 

dt 


dA  ,  dA  , 
sd  ed 


dt 


dt. 


(41) 


From  the  geometry  we  get 


A  =  A  -  cos 
sd  pa 


-1 


cos  (v  -  V  ,)/c03  X 
s  sd  • 


(42) 


and 


dA 


sd 


dt 


cos  (v  -  V  ,)  sin  X 

3  sd _ pfi 


pa 


sin  (a 


pa 


A  )  cos  X  ^ 
sd'  pa 


(43) 


Again  using  small -angle  approximations,  as  in  the  derivation  of  Eq,  (9), 
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and  Eqs.  (3$)  and  (36),  ve  substitute  Into  Eq.  (43)  and  get 


dA 


sd 


(U4) 


Hie  second  term  on  the  right  aids  of  Eq.  (4U)  is  amaii  conqpared  to  one 

except  for  values  of  (a  -  A  )  near  180°,  since  the  product  X  1  is 

pa  sd  pa  p 

small  for  all  planets  of  interest.  For  (A  -  A  .)  I  180°  the  inclination 

pa  sa 

of  the  transfer  orbit  plane  la  approximately  90°.  transfer  orbits  of  this 
inclination  would  require  extremely  high  velocities;  therefore  this  term 
is  neglected  and 


Hie  time  rate  of  (^nge  of  the  Earth's  heliocentric  longitude  is  given 


(45) 

by 


Substitution  of  Eqs.  (4^)  and  (46)  into  Eq.  (41)  gives 


(46) 


(47) 
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D5RIVATI0N  OF  (U)  OF  METOOD  B.3 


d(Atg)  d[tg  -  tg  (fixed)]  dt^ 


^*d 

Since  tg  =  f(ng,  n,  n^) 

dt  dt  dn 

_ s  s  s 

dt ,  “  Sir  dt - 
d  s  d 

dt  , 
dq  dt^ 

bt 

* 

dll 

(48) 

'nie  total  derivatives  are 

as  follows: 

"d 

3 


(49) 


If  we  assume  that  r  ,  and  r _  are  constants  and  use  the  fact  that  v_ ,  is 

sd  pa  sd 

fixed,  then  using  Eq.  (l) 


da 

_ E 

dt. 


^s  ^ 


da  de  dv 

S  S  8 

de  dv  dt. 

s  s  d 


(50) 


since  e^  -  f(r^^,  v^,  v^).  Using  Eq.  (3O) 


dv  dv^  "^Tja  ^s  ^na 

“i  ■  ^  "d  *  ^  “d 


(51) 


dv  dX  dv  dA 

8  _]£  ^  8  _JB 

dX  dt  dA  dt 

pa  a  pa  a 
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since  t  =  t  +  t  emd  t  is  constant, 
ads  s 

ihe  partial  d£;rivative  ov  /hx  contains  the  sin  X  ,  and  the  total 
^  3  pa 

derivative  dA,  /dt  contains  tan  i  .  Since  both  X  and  i  are  sdbJJ.,  the 
pa'  a  p  pa  p  ^ 

product  of  these  derivatives  is  neglected  in  Eq.  (5I)  and  subsequent  steps 
of  the  derivation.  Thus,  Eq.  (5I)  becomes 


dv  c)v  dA 

_ 5  _  a  pa 


dt 


OA  dt 
pa  a 


(52) 


and  Eq.  (4-9)  becomes 


dt. 


-a  / 

V 


:  \q  C  rt\r 

S  '  '[S  '~y  S  • 

^5/2  Se  'ov  dA  ‘'pa 

2a  '  3  s  pa 

8 


(53) 


According  to  the  assumptions  above  and  Eq.  (19) 


dn  ^ 

dt ,  ~  8a 


da 


dt 
s  a 


dt 


a 


da  X  dv 

s  ,  on  oc  s 

dT  ^c  ST  dT 
3  a  E  a 


(54) 


Substitution  of  Eqs.  (50)  and  (52)  into  Eq.  (54)  gives 


dt . 


da  c3e 

+ 

da  de  ov 
3  s  s 


8c_\  ^^s  • 

Sv  j  8a  pa 
sy  pa 


(55) 


By  using  Eq.  (20)  in  a  similar  manner 


dt ,  \  da  8e  8v 

d  \  3  G  0 


^1  Sc  \  • 

Tk"  ST'  S”  V 

a  /  pa  ^ 


(56) 
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By  using  Eqs.  (I9),  (20),  and  (2I)  wc  find  that 


L  tan  V2  ^  = 

s 


1  .  In 

XT  =  -  r 


c  2a  sin  n 
s  ' 


2a  _  sin 


(56a) 


=  -  r  r  sin  (v  -  v  .) 
ov^  c  sd  pa  s  sd 


By  solving  Eq.  (l)  for  e  and  then  for  a  we  get 

6  S 


r  a  -  1' 
sd  pa 


and 


r  cos  V  -  r  ,  cos  v  , 
pa  s  sd  sd 


"-Bd  \  td' 


1  -  e‘ 


The  partial  derivatives  are 


Oa 

8 

8 


a  (1  +  e  )  -  r  , 
s  _ s  sd 

e  (1  -  e-) 

3  3 


and 


c)e  e”^  r  sin  v 

3  8  pa  s 

^  ""r 

s  sd  pa 
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Using  Eq.  (30)  we  get 


substituting  for  the  partial  derivatives  of  Eqs.  (53)»  (55)» 

(56)  and  then  substituting  these  equations  into  Eq,  (U8)  and  simplifying, 
ve  get 


If  ve  eliminate  the  remaining  partial  derivatives  using  Eqs.  (25)»  (26),  and 
(27)  and  simplify,  we  get 


S  t  a  .  cos  n)  tan  n/2 

'pa\V  (1-^  (^d-  V  /'■ 


t  (1  =  cos  tan  n^/2]  ♦  ^  Bin  (v^  -  v^^)  i  tam/2  ^  tan  r^j2 


(57) 
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Tlie  multiple  signs  in  bre.c.kets  depend  on  the  transfer  orbit  as  foiJ.O'V're: 

Direct  =  -,  +  and  +,  +  Perihelion  =  -  and  +,  - 

Aplielion  =  +,  +  and  +  Indirect  =  +,  -  and  - 

DERIVATION _0Ii_E-vi_.  (12)  of  METHOD  C.l 


dt 

a 


d  a  (min) 
s  _ _ 

dt 

a 


(58) 


since  for  fixed  V  ,  a  choice  of  t,  determines  r  ,  and  consequently  a  . 

sd  d  ed  ^  ^  s 

Thus 


da 


s 


dt 

a 


=  0 


The  equation  for  a^(min)  is 


a  (min) 
s 


(r 


sd 


+  r 


pa 


c) 


(59) 


where  r  ,  =  r  We  assume  that  r  is  constant.  As  a  result  of  this 
sd  ed  pa 


da  (min) 
s 

dt 

a 


ha  (min)  , 

s)  '  dc 

5c  dt 

a 
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Using  Eq.  (2)  for  c  and  Eq.  (3O)  we  get 


Since  the  coefficient  of  contains  the  sin  and  X^^  contains  tan  i^, 
the  second  term  in  the  bracket  is  neglected  as  in  the  derivation  of  Eq.  (9). 
Equation  (58)  may  now  be  written  as 


da^(min) 

dt 

a 


(61) 


^  using  Eqs.  (21),  (3O),  and  (59),  the  partial  derivatives  are  evaluated 
and  then  substituted  into  Eq.  (61).  After  simplification  we  have 
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1316  derivative  dt  (step  4}/dt  =  1  since  for  t, (fixed)  a  change  in  the 

S  £l  Q 

chosen  value  of  t^  causes  an  equal  change  in  t^. 


dt  ^t  dn  8t  ,  f^t  dn^ 

_ 3  _  a  _ s  ^  s  di)  3  _1 

dt  dn  dt  Sri'  dt  SiiT  dt 

n  o  O  •  Q  *  I 


1  a 


For  fixed  V  ,  and  t.,  a  is  determined  and 
sd  d^  s 


3/2 

dn^  d(  yiT/a^  ) 


Using  Eqs.  (19),  (2O),  (2I),  and  (60)  we  get 


dr)  _  ^  dc  _  ^  8c 

dt  "  Sc  dt  "Sc  8(v  -  V 


3  ad 


8(v  -  V  ,) 
8  ad 


^  (64) 

pa 


^  dc_  Sc  ^‘''a  -  •  ,,,, 

dt  **  Sc”  dt  Sc  S(v  -  V  ,)  Sa  pa  ^ 

a  a  '  B  ad'  pa 


lity  replacing  tlie  total  derivatives  and  simplifying,  Eq.  (62)  becomes 


At  )  -s  S(v  -  V  ,)  .  X,  St  St),  St 

_ ^  ^  Sc  ^  s  sd'  ^  ^  a  ^  '1  3  ^ 

dt  °  S(v  -  V  , )  Sa  'pa  Sc  Sn  Sc  SnT 
a  8  sd'  pa  ^  '  '1 


<i(Aty 


r  n  ./p  P 
sd  V  s  r 


St  Sr)  St  ' 

"  ^dn  S  Sir  ST  Sfi^ 


~  -  1  (66) 
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The  partial  derivatives  of  Eq.  (66),  \rtiich  are  given  by  Eqs.  (26),  (27), 
and  (56a),  are  replaced,  and  we  get 

d(Atg)  r^^  ,  + 

dt  “  2c  r  ^  ■  '^sd^  -  tan  n^/2)  -  1  (67) 

The  signs  of  the  terms  in  parentheses  depend  on  the  type  of  transfer  orbit 
as  follows: 


Direct  =  +,  + 


Perthellon  =  +,  - 


Aphelion  =»  -,  + 


Indirect 


DERIVATION  OF  E.^.  (16)  OF  METHOD  D.l 


^y  choosing  t^  in  addition  t/:  ‘■.he  parameter  t^,  the  t^  is  determined, 

and  consequently  (v^  -  v^^)  is  a  constant.  Again  we  assume  that  r  is 

P® 

constant. 


dtg(step  9) 


^  ^  aa_  .  iij. 


The  total  derivatives  are  written  as  follows: 


372 


(66) 
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8c  &  r  sin  v 

s  s  pa  s 

^  =  -  r‘  -  r' 

s  sd  pa 


e  r  ,  3in  v  , 
s  sd  _  sd 

r  ,  -  r 
3d  pa 


da  e  r  ,  sin  v  , 

_ ^  B  sd  sd 

8v  “  2 

sd  1  -  e 

3 


Eq.  (70)  becomes 


da  e  a(l  +  e^)  -  r  , 

- -  = - - - sin  V  -  r  ,  sin  v 

dv  ,  ,  2  r  ,  -  r  pa  s  sd  sd' 

sd  1  -  e  sd  pa 


r  ,  sin  V  . 
sd  sd 


The  partial  derviatives  of  Eq.  (69)  are  given  by  Eqs.  (23)  -  (27)- 
Using  these  equations  and  Eq.  (7l)>  Eq.  (69)  can  be  written  as 


dt  / a  e  a  (1  +  e^)  -  r  , 

s/ss  s'  s'sd/  ,  ,  \ 

_ —  _  /  —  - - y  -  V  ,  sin  V  -  r  ,  sin  v 

dv  /  p  2  r  ,  -  r  pa  s  sd  sd  sd 

sd  V  3  1  -  e  sd  pa 


^  +'“nt  -(l~  cos  q)  ton  q/2  -  (l  -  cos  q,  )  tan  q,/2 

Cl  S  S  '  * 
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The  multiple  signs  in  the  second  brackets  depend  on  the  type  of  transfer 
orbit  as  follovs: 


Direct  =  -,  + 

Perihelion 

Aphelion  a  +,  + 

Indirect  = 

RM-2881 
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